Derivation of Static Low-Energy Effective Models by ab initio Downfolding Method 
without Double Counting of Coulomb Correlations: 
Application to SrV0 3 , FeSe and FeTe 
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Derivation of low-energy effective models by a partial trace summation of the electronic degrees of 
freedom far away from the Fermi level, called downfolding, is reexamined. We propose an improved 
formalism free from the double-counting of electron correlation in the low-energy degrees of freedom. 
In this approach, the exchange-correlation energy in the local density approximation (LDA) is 
replaced with the constrained self-energy corrections defined by the sum of the contribution from 
eliminated high-energy degrees of freedom, E^ and that from the frequency-dependent part of the 
partially screened interaction, AE L . We apply the formalism to SrV03 as well as to two iron- 
based superconductors, FeSe and FeTe. The resultant bandwidths of the effective models are nearly 
the same as those of the previous downfolding formalism because of striking cancellations of T, H 
and AE L . In SrVC>3, the resultant bandwidth of the effective low-energy model is 2.56 eV (in 
comparison to LDA: 2.58 eV, full GW approximation: 2.19 eV). However, in the non-degenerate 
multi-band materials such as FeSe and FeTe, the momentum dependent self-energy effects yield 
substantial modifications of the band structures and relative shifts of orbital-energy levels of the 
effective models. The FeSe model indicates a substantial downward shift of the x 2 — y 2 orbital 
level, which may lead to an increase in the filling of the x 2 — y 2 orbital above half filling. It 
suggests a destruction of the orbital selective Mott phase of the x 2 — y 2 orbital in agreement with 
the experimental absence of the antiferromagnetic phase. 
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PACS numbers: 71.15.-m, 71.28.+d, 71.10.Fd 



I. INTRODUCTION 



Strongly correlated electron systems have been a cen- 
tral subject of condensed matter research for decades. 
A challenge is to unveil their electronic structures from 
first principles. A major advance has recently been made 
in this direction^ by taking advantage of the hierarchi- 
cal energy structure of strongly-correlated-electron ma- 
terials, which commonly have only a small number of 
bands near the Fermi level in contrast to usual highly 
dense band entanglement at energies far away from the 
Fermi level. In this general framework of multiscale ab 
initio scheme for correlated electrons (MACE) 1 -, global 
dense and entangled band structure is calculated based 
on density functional theory (DFT) 3 -^. The degrees of 
freedom far away from the Fermi level is traced out per- 
turbatively in the spirit of the renormalization group, 
which generates an effective low-energy model consist- 
ing of a small number of bands near the Fermi level, 
which we call "target bands" . This downfolding pro- 
cedure is efficiently performed by the constrained ran- 
dom phase approximation (cRPA)^^ after constructions 
of the maximally localized Wannier orbitalsi^. Derived 
ab initio low-energy models are represented in general 
by multi-band fermions on lattices, which are typically 
extended-Hubbard-type models. The low-energy mod- 
els arc solved by high-accuracy solvers such as vari- 
ational Monte Carlo (VMC) 53 , path-integral renormal- 



ization group (PIRG) 10 ' 11 , dynamical mean field the- 
ory (DMFT^ii 3 - with its cluster extensions 2 * 1 ^ - — and 
functional renormalization group (FRG)^. The validity 
of the cRPA is ascribed to the small vertex correction 
thanks to the above hierarchical structure 1 -, while an ex- 
tension of the cRPA by using the FRG has also been 
proposed 2 ^. 

This hybrid scheme referred inclusively as MACE 
has been applied to a wide variety of materials; 
semiconductors 2 ^, transition metals&2Zr— , transition- 
metal-oxide o 28 ' 30 " — , transition-metal-oxide interface 4 ^, 
molecular organic conductors 4 ^—, and Fe-based layered 
supcrconductor o 45 ' 46 by using the solvers with the choices 
of VMC 4 !^, DMFTia-^i and the Ab initio ef- 

fective models were derived for other challenging complex 
materials with a large number of atoms on a unit cell as 
weU&S&. 

Quantitative accuracies of obtained physical quanti- 
ties in comparison to experimental results support the 
validity of the scheme. However, there remains a ba- 
sic question about the counting of the electron correla- 
tion: The correlation effects are, though insufficiently, 
taken into account, in principle, in the DFT calculation 
as the exchange correlation energy. Since the correlation 
effects contained in the effective low-energy model are 
considered when the model is solved by the low-energy 
solver, there exists a well-known double counting prob- 
lem. Removing the Hartree contribution can be consid- 
ered easily when solved by the low-energy solver—. On 
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the other hand, the double counting arising from the ex- 
change correlation has simply been assumed to be small 
in most of the studies in the literature. Removing the 
double counting as well as removing the influence of the 
original LDA/DFT bias is a desired step for the consis- 
tent and ab initio descriptions of the correlation effects 
in the present downfolding scheme. Several attempt in 
this direction has been put forward in the literature^—. 
Since it is not a practical choice to completely throw away 
the DFT/LDA calculations in obtaining the global band 
structure, it is better to remove the bias and the double 
counting as a correction as much as possible within the 
present framework. 

If the correlation effects are taken into account par- 
tially in the exchange correlation in the LDA calculation, 
narrowing of the bands coming from the effects within the 
target bands, however small it is, is regarded as the dou- 
ble counting, and the resultant parameters underestimate 
the transfer integral in the effective model. 

To overcome the double-counting problem, in the first 
step of this paper, we subtract the whole exchange cor- 
relation energy in the LDA to remove the double count- 
ing and to completely eliminate the bias of the LDA. 
Instead of the exchange-correlation potential V xc , we 
then add a constrained self-energy arising from the high- 
energy region Y, H = —G H W based on the GW approx- 
imation (GWA) 63 i 64 by excluding the self-energy contri- 
bution from the low-energy target degrees of freedom. 
Here, G H is the Green's function for the degrees of free- 
dom outside the target bands (high-energy bands) and W 
is the fully screened interaction in the RPA as we define 
below. 

In general, the one-body parameter calculated from 
such constrained GWA and the two-body parameter cal- 
culated from the cRPA have frequency dependences orig- 
inating from the constrained self-energy H H {w) and the 
constrained polarization P H (ui), respectively. Although 
the frequency dependence of the interaction has been con- 
sidered in the literatur o 39 i 54 ' 65 , in practice, it is difficult 
to treat the dynamical Hubbard model by keeping both 
the frequency and momentum dependences of fluctua- 
tions at high accuracy. So far, although DMFT is able to 
incorporate frequency dependences of the effective inter- 
action, it is actually difficult to consider spatial fluctua- 
tions accurately in DMFT. Therefore, in this paper we 
make an attempt to derive reliable static effective fermion 
models on lattices. 

Then, in the second step of the present study, we renor- 
malize the frequency dependence of the effective interac- 
tion to derive static low-energy effective models. This 
Hamiltonian formalism allows the consideration of spa- 
tial fluctuation effects accurately by low-energy solvers. 
We renormalize the frequency dependence of the con- 
strained self-energy T, H into static one-body parameters 
by using the renormalization factor Z H , namely by using 
the linear frequency dependence of the real part of T, H . 
As we will see later, the constrained self-energy has a 
well-behaved and gentle frequency dependence near the 



Fermi level, and thus the expansion around the energy 
eigenvalue is justified. In addition, we renormalize the 
effect of frequency dependence of the partial screened 
Coulomb interaction W p (u>) into the one-body part as 
the self-energy^. 

We next apply the present formalism to several typ- 
ical examples: First we apply to SrV03, where the ti g 
orbital degeneracy is retained. Second we apply to FeSe 
and FeTc, typical iron-based superconductors, where the 
lift of the orbital degeneracy introduces additional com- 
plexity. 

In these examples, we show that after removing the 
double counting in the low-energy space obtained in the 
first step, the bandwidth and the transfer integral in- 
crease typically about 30-50% than those of the LDA. 
When we perform the second step and renormalize the ef- 
fect of the frequency dependence of the effective Coulomb 
interaction screened by the high-energy bands into the 
mass enhancement, then the bandwidth and the transfer 
integral remarkably recovers to those of the LDA. 

However, in the case of the multi-band systems with 
non-degenerate orbitals, the removal of the double count- 
ing causes the orbital dependent shift of the chemical 
potential. Because the constrained self-energy is not 
only frequency dependent but also wavenumber and or- 
bital dependent, the low-energy band after excluding the 
double-counting is different from that of the LDA, es- 
pecially in non-equivalent multi-orbital systems such as 
Fe-based layered superconductors. However, A warning 
is that this band structure should not be compared with 
experiments, because the experimental band structure is 
a consequence of the many-body effects within the tar- 
get band and the direct comparison is possible only after 
solving the low-energy effective model. 

In Sec. 2 we propose the formulation. Section 3 de- 
scribes the global band structures together with the 
derived effective models after the improved downfold- 
ing for the examples of a transition-metal-oxide SrV03 
(SVO) and the Fe-based layered superconductors, FeSe 
and FeTe. Section 4 is devoted to discussions and sum- 
mary. 



II. METHOD 

A. Constrained Random Phase Approximation 

In the LDA, the band structure is calculated from the 
Kohn-Sham equation 

w LDA^DA (r) 

= [-Iv 2 + V cxt (r) + V H (r) + y xc (r)]^ LDA (r) 

= e ^ D >), (I) 

where V cxt is the potential of atomic nucleus, V K is the 
Hartree term of the Coulomb interaction, and V xc is the 
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exchange-correlation potential. The energy eigenvalue is 
denoted by e^ DA , and ipf^ is its eigenstate. 

Now we derive the effective model for the target bands 
as a typical case. Our target bands here are the d elec- 
trons. From the eigenvalues and eigenstates of Eq. (JTJ) , 
the transfer integral of the low-energy effective model is 
given by 



,LDA 



(R) 



n 



LDA 



(2) 



where 4>mR is the maximally localized Wannier function 
(MLWF) of the mth orbital localized at the unit cell R in 
the low-energy space derived from a linear combination 
of the Kohn-Sham (KS) wave functions by following the 
conventional prescription^. 

The partially screened Coulomb interaction for the d 
bands is given by 



U mn (R,cj) = 



W p 



1 - V P H ' 



(3) 



(4) 



where v is the bare Coulomb interaction, and P H is the 
polarization without low-energy L-L transition, P H = 
P — P L . Namely, P H is defined from the total polar- 
ization subtracted from P L , the polarization contributed 
only from the L space. The indices "L" and "iZ" express 
the low- (within target bands) and high-energy (including 
outside of the target bands) parts of the Hilbert spaces, 
rcspectivelyii The frequency dependence of W p comes 
from the frequency dependence of the high-energy polar- 
ization P . The partially screened Coulomb interaction 
is usually written as W r in the literature, but, in this 
paper, we express it as W p to avoid confusions with the 
renormalized Green's function in the low-energy space. 



B. Double-Counting of Self-Energy in Low-Energy 
Space 

The transfer integral in Eq. @ formally contains the 
self-energy of the low-energy space in the form of the 
exchange-correlation energy V xc in the DFT/LDA. How- 
ever, the parameters of the low-energy effective model 
should not contain such an energy of the low-energy 
space, because the self-energy effect in the low-energy 
space should be considered when one solves the low- 
energy effective models. Therefore, this contribution 
has to be subtracted, but it is not easy to separate 
V xc (n(r) a ii) to the high- and low-energy parts. 

On the other hand, in the GWA beyond the LDA, we 
can easily divide the whole self-energy GW into the low- 
and high-energy parts: 



GW = G L W + G H W, 



(5) 



where the whole Green's function G is given by the sum 
of the low- and high-energy propagators, G L and G H , 



respectively, and W is the fully screened Coulomb inter- 
action. 

When one subtracts the exchange-correlation energy 
V xc from the one-body energy in the Kohn-Sham LDA 
Hamiltonian H LDA and add a constrained self-energy 
Yi H = —G H W to that, we can calculate the transfer in- 
tegral by excluding the self-energy contributed from the 
low-energy space only: 



\H WA -V™ + V«\rt R ). (6) 



This allows to eliminate the double counting of the self- 
energy in the low-energy space when one solves the low- 
energy model honestly. 

Moreover, it allows a formalism to calculate the one- 
body parameter beyond the DFT/LDA schemes. The 
good convergence of the GW scheme in terms of the 
perturbative expansion and the irrelevance of the vertex 
correction^ assure the accuracy of 



C. Static Hubbard Model 



In the previous section, we consider the "dynamical" 
effective Hubbard model 



% cff — t^ na (Ri - Rj,tj)d ina dj m <T 

ij rana 

+ g ^ ] ^ ] i U mn<J p(Ri — Rj , uj)dl niJ dj mp dj m pdi na 

ij mnap 

~\~ Jmnapi^Ri Rj : ^) ^d^ na dj Tn pdinpdjrac7 

+ di na d inp dj mp dj m cr) >, (7) 



where d\ na idina) is a creation (annihilation) operator of 
an electron with spin a in the nth MLWF centered at Ri . 
The effective Coulomb interaction U and the exchange 
interaction J are given by 



\w p \^ R ^ R 



and 



J mn (R) = (^„4l^lMii> 



L iL 



(8) 



(9) 



respectively. 

The renormalized Green's function of the low-energy 
space G Lp , where the constrained self-energy from the 
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high-energy space is implemented, is written as 



G Lp (tu) = 



G L0 



i - g lo i: h 
1 



(#LDA _ y x 



E*( W )) 
1 



(10) 
(11) 



where the static one-body part t H is given by 

tmni R ) = {<t>mQ\H H \<l>nR), 



(tf LDA - V xc + S ff (e LDA ) + & w=eLDA (u; - e LDA )) 



(12) 



Z ff (e LDA ) 



oj - (i? LDA + Z«(e LDA )(-y xc + S ff (e LDA ))) ' 

where the "bare" Green's function G L0 is 

1 



(13) 



G lo (uj) = 



uj - (H LDA - V xc ) ' 
and Z ff (e LDA ) is the renormalization factor of T, H : 

<9Re£ ff 



Z H (e) = 1 - 



eta 



(14) 



(15) 



We note that 7J LDA is not renormalized by because 
the Fermi level must be fixed by the electron number 
and we have taken into account the uj dependence around 
e LDA . Because T, H = —G H W does not contain the poles 
of G L , Z H is close to 1 in contrast to the substantially 
small renormalization factor obtained for the full self- 
energy GW, as we will see later. 

If we can ignore the imaginary part of T> H and Z H in 
the numerator of G Lp (u>), and the terms higher than the 
linear one in the u> dependence in the real part of T, H , 
we obtain 



G Lp (uj) 



UJ 



(H LDA + Z H (e LDA )(^V^ + Re£ ff (e LDA ))) ; 



(16) 



or, in the case where the self-energy is expanded around 
the Fermi level instead of the eigen value in the LDA, 



G Lp ( 



1 



- Z H (0)(H LDA - V™ + Re£ ff (0)) 



, (17) 



which preserves the number of the electron in eq. (|13p 
and has the form of the "dynamical" effective Hubbard 
model, containing the uj dependence in the two-body 
part, 



J- j ma 



ij mricr 



ij mncrp 



Jmnap (R Rj , k>) (d\ na dj mp di n pdj m(7 



+ 4na4np d jmpdj m a)\, (18) 



(19) 



n H = n LDA + z n (e 



-ff /,LDA 



= Z H (0)(H 



)(-V xc + RcT, H (e hDA )) 

Re£ ff (0)). (20) 



LDA 



- V s 



In Eq. ((20)) , we employed the renormalization factor Z H 
only with multiplying by the diagonal part of RcS ff , 
whereas Z H for the off-diagonal elements is fixed to unity. 
Eigenstates and % H are determined self-consistently. 

Actually, the u dependence requires the Lagrangian 
description instead of the Hamiltonian description. How- 
ever, if the frequency dependence is small in the energy 
range of physics to be considered by the low-energy effec- 
tive model, one may ignore the frequency dependence by 
replacing U (uj) and J(u>) with its static limit ui — > 0. This 
is indeed the case when the low-energy bands are well sep- 
arated from the high-energy ones as is normally expected 
in the strongly correlated electron systems^. Then one 
may employ a static value 



W p (uj = 0) 



1 - vP H (0) 



(21) 



instead of Eq. (gj) . 

To further renormalize the effects of w-dependence of 
P H in the GWA, we append the self-energy correction. 
A choice of this correction is 



AS 



HF 



G Lp (W p (uj) - W p (0)). 



(22) 



Note that W p (uj) ~ W p (0) vanishes at w = 0. Equa- 
tion (|22|) is the "Hartree-Fock" approximation (HFA), 
because G Lp is the "bare" Green's function in the sense 
of the absence of the self-energy arising from L-electron's 
many-body effects and W p is the "bare" Coulomb inter- 
action in terms of the low-energy effective model. 

In the case of the GWA beyond the HFA, the self- 
energy correction is calculated as 



A^ W = G L (W -W u ), 



where 



W 



W p (0) 



1 - W p (0)P l (uj) 



(23) 



(24) 



Here W u (uj) is equal to the fully screened interaction 
W(u>) at oj = 0, while W u (uj) approaches W p (0) for ui — > 
oo. On the other hand, W approaches v in the ui — > oo 
limit. In this one-shot GWA, the self-energy expected 
in the level of the GWA of the effective model, namely 
G L W U , has been subtracted from the full GW self-energy 
G L W, since G L W U should be considered when one solves 
the model. In this paper, we calculate AS L by the one- 
shot GWA according to Eq. ([2"3"|). 
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Thus, the renormalized Green's function for the static 
Hubbard model 



G LpA (cu) = 



G Lp 



1 - G l pAE l 



1 



ui - (H LDA - V* c + Ti H 



AS L ) 



(25) 
(26) 



- {H LDA + Z HA {e LDA )(-V* 



<9(Rc£ ff + ReAS L ) 



(28) 



is obtained. 

Similarly to the reduction to Eq. ()16[) . even with 
this AE L correction as well, we reduce the renormalized 
Green's function with Z factor to 

G LpA {co) = 

1 

- _ (#LDA + Z ffA( £ LDA)(_y X c + Rc (vjff + AE i )(e LDA ))) 

1 

~ uj - Z HA (0)(H hDA - V* c + Re(£ ff + A£ L )(0)) ' 

(29) 

In the practical calculation below, we employ AS^J F in 
Eq. d22J) for AS L . Then the "static" effective Hubbard 
model, where all the effects of frequency dependences in 
the one- and two-body parts are incorporated into the 
renormalization of the one-body part as 



M jrna 



ij mna 
ji m.ri.rtn 



tj mncrp 



Jmnap{Ri Rj > fy{d\ na d^ mp di np dj 



+ 4na4np d jmpdjma) \ , (30) 



where the renormalized static one-body part t HA is given 
by 



(31) 



,LDA + Z HA (£ LDA )( _ y xc + Rc ^H + A E i )(e LDA )) 

= Z HA (0)(n LUA - V xc + Re(S H + AE L )(0)). (32) 



The renormalized effective models of Eqs. ([20]) and 
(|52"|) with Eqs. (HU) and (JSU still contain a double 
counting of the Hartree term arising from the low-energy 
(target) degrees of freedom in the one-body part. This 



TABLE I: Band width of the i2 9 band in each calculation of 
the self-energy. Units are given in eV. The band width of 
"H H " is calculated according to Eq. |16J) and that of U HA 
is calculated according to Eq. (|29p respectively. The band 
width except for "LDA" includes the tj-dependence effect of 
the self-energy through the renormalization factor. 



SVO 



(Tj h + AE 1 )(e LDA ))) ' 
(27) 



width 



LDA GWA H H H 



2.58 2.19 3.41 2.56 



Hartree double counting may easily be subtracted in the 
model^, because the Hartree term in the LDA and that 
in the model are essentially the same. For example, to 
subtract the double counting of the Hartree term for the 
on-site potential, one should employ the one-body part 



.dcf 



.{Ri) — t" A a {Ri) 



E 

jnp 



u„ 



){Ri 



Rj,0){n L 



jnpl 



where the average (n^ np ) 



(33) 

is calculated from the Hartree 
approximation of H LDA combined with the partially 
screened interaction calculated from cRPA as listed in 
Ref4£. This prescription is based on the observation that 
the Hartree potential originated from the low-energy tar- 
get degrees of freedom considered in the LDA may well 
be reproduced in the Hartree approximation of the pre- 
vious model obtained without considering the self-energy 
effect, because the Hartree potential must be essentially 
the same. Then we can determine the on-site potential 
tmmcr(Ri) lree & om the double counting. 



D. Computational Conditions 

Computational conditions are as follows. We calculate 
the band structures of the transition-metal-oxidc SVO, 
and the Fe-based layered superconductors FeSe and FcTe 
based on the DFT/LDA^. The band structure calcula- 
tion is based on the full-potential linear muffin-tin or- 
bitals (FP-LMTO) implementation^. The cRPA and 
GW calculations use a mixed basis consisting of products 
of two atomic orbitals and interstitial plane waves^. In 
the LDA calculation, 8x8x8 fc-mesh is employed for 
the SVO, 12 x 12 x 6 fc-mesh is employed for FeSe and 
FeTe. In the cRPA and GW calculation, 6x6x6 fc-mesh 
is employed for the SVO, 3x3x3 fc-mesh is employed 
for FeSe and FeTe. 



III. RESULT 
A. SrVOa 

In this section, we derive the low-energy effective model 
for the V ti g bands of the transition-metal-oxide SVO. 
Figure [T] (a) shows the band structure of SVO in the 
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(a) SVO LDA 
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FIG. 1: (Color online) Electronic band structures of SVO 
obtained by (a) the LDA and (b) the GWA. The zero energy 
corresponds to the Fermi level. In the GWA, the self-energy 
is calculated only for the t2 g Wannier band. 



V h g xy 




FIG. 2: (Color online) Isosurface of the maximally localized 
Wannier function for ±0.05 a.u. for the V d xy orbitals in 
SVO. 



DFT/LDA. SVO is a paramagnetic metal^ and has a cu- 
bic structure with lattice parameter asvo = 3.843 ASS. 
The three conduction bands at the Fermi level are de- 
rived from the ti g orbitals of V sites, where the octa- 
hedral crystal field of O -2 partially breaks the five-fold 
symmetry of the 3c? orbitals into the lower orbitals of the 
ti g and the higher orbitals of the e g . The band width 
of the t2 9 bands is 2.58 eV. One conduction electron per 
unit cell is accommodated in the ti g bands. Hereafter, 
we regard such three conduction bands of the ti g orbitals 
as the low-energy bands and the rest of the whole band 
structure of SVO as the high-energy bands. The top of 
the occupied high-energy band is the O 2p band and the 



bottom of the unoccupied high-energy band is the V e g 
band. 

Figure Q] (b) shows the band structure of SVO in the 
GWA, where we calculate the self-energy only for the V 
tig orbitals after removing the exchange correlation en- 
ergy in the LDA. Figure [D shows the isosurface of the 
Wannier function of the V d xy orbital, where the V d xy 
atomic orbital hybridizes the nearest O p atomic orbitals. 
The Fermi energy of the low-energy bands in Fig.[T](b) is 
determined by the occupation number in the low-energy 
space separated from the whole Hilbcrt space. Due to 
the large self-energy mainly from the low-energy space of 
the tig bands, the band width in the GWA is about 15% 
smaller than that in the LDA (see Table H]). In Fig. |31 
the cj-dependence of the self-energy of SVO in the one- 
shot GWA is not smooth near the eigenvalues of the V 
tig bands. Incidentally, when one chooses [—10, 10] (eV) 
instead of [—0.5, 0.5] (eV) as the u range to calculate the 
renormalization factor, then the average of Z consider- 
ably increases from 0.55 to 0.74. This is an indication 
that the linear approximation of the self-energy collapses 
in the low-energy band near the Fermi level. 

To calculate the transfer integral without double- 
counting terms in the low-energy space, we calculate the 
constrained self-energy originated only from the high- 
energy space, Y. H = —G H W, instead of the exchange 
correlation term y xc in the LDA or the full self-energy 
GW in the GWA. We show the w-dependence of the full 
self-energy GW and those of the constrained self-energies 
G L W and G H W in Fig. [3], where the real and imaginary 
parts of the correlation part at T-point are presented in 
Fig. [3] (a) and (b), respectively. As one can see in Fig. |3l 
the w-dependence of GW mainly originates from G L W, 
and thus the w-dependence of G H W is very smooth and 
weak compared to those of GW and G L W and the par- 
tial renormalization factor Z H is close to 1 (see Table HI)) . 
These results arc obtained from the lst-shot GWA, the 
first-step of the iterative approximation of the GWA from 
the DFT/LDA. The smoothness of G H W originates par- 
tially from the validity of G H as the initial-state of the 
iterative calculation of the GWA. Because of the smooth- 
ness and the small frequency dependence of G H W, the 
truncation up to the first-order expansion of Yj H in fre- 
quency in Eq. (fT2")) and the neglect of ImS H and Z H in 
the numerator of G Lp in Eq. (fTB|) arc justified. Actually, 
Z H in this calculation does not sensitively depend on the 
choice of the w range. Figurc[4]is the w-dependence of the 
transfer integral without double-counting terms to the 
nearest [1,0,0] and next-nearest [1,1,0] V ti g orbitals. 
Similar to the self-energy of the ti g Bloch function in 
Fig. [3j the w-dependence of the transfer integral without 
double-counting terms is very weak. In the next para- 
graph, the effect of the weak w-dcpcndcnce is renormal- 
ized to the transfer integral through the renormalization 
factor Z H . As one can sec m the list for Z H in Tabic M 
the band narrowing effect by the u> dependence of Y, H is 
small and nearly uniform in fc-space. 

Figure [5] shows the band structure after considering 
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FIG. 3: (Color online) (a) Real parts and (b) imaginary parts of correlation self-energy of the t^g Bloch functions at F = (0, 0, 0) 
in the Kohn-Sham gauge. The zero energy corresponds to the Fermi level. Inner window is the magnification only for "G^W". 
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FIG. 4: (Color online) Dynamical transfer integral of the tig orbitals of the V sites without double-counting terms: t (cj) = 
((^1% - F xc + E \<j> L ) (see Eq. @). The left panel is the nearest hopping [1, 0, 0] of the V t 2g orbitals and the right panel 
is the next-nearest hopping [1, 1, 0]. These values are dynamical and, therefore, do not include the effect of the renormalization 
factor Z H . 



the self-energy G H W, where we calculate G H W for V 
t2 g bands instead of V KC and subtract the constant shift 
of the self-energy from the low-energy bands. The band 
structure of the t2 g states is calculated according to Eq. 
([TTJ)) . where the effect of the renormalization factor Z H 
is included. In Fig. [5j the Fermi level for the low-energy 
band is determined from the occupation number in the 
low-energy space. On the other hand, because the energy 
is determined not only from the low-energy space but also 
from the high-energy space, the energy is calculated from 
all dcgrcc-of-frecdom including the high-energy space. If 
we calculate the self-energy from the low-energy Green's 
function G L in the GWA (G L W), then the band struc- 
ture in the low-energy space corresponds to that in the 
full GWA in Fig. ED (b). 

We summarize the values of the transfer integral in Ta- 



ble Hill The transfer t H in the left column of Table Hill 
corresponds to Fig. [51 where the w-dependence of Y, H is 
renormalized to Z H . The renormalized values obtained 
further from the cj-dependence of the partial screened 
Coulomb interaction W H (oS), namely, AS L and t HA , are 
discussed in the next paragraph. We note that the pa- 
rameter of the t2 g orbital has the same symmetry with 
the V t2 g orbitals in the real SVO. In the LDA, the largest 
value of the nearest hopping is between the same symme- 
tries t^P^ y (l, 0, 0) and is about 270 meV, whereas that of 
the next nearest hopping t^^ y (l, 1,0) is about 90 meV. 
Compared to the results of the LDA and the full GWA 
in Fig. [I] the band width of the V t2 g in Fig. [5] is much 
larger because of the absence of the self-energy originat- 
ing from the low-energy space G L W (see Table [T|. Corre- 
spondingly, the transfer integral for the nearest-neighbor 
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TABLE II: Renormalization factor Z for each V t^ g Wannier band in bulk SVO. Z is the partial renormalization factor: Z — 
(1 — d g |^ =0 ) -1 ■ "Average" is the average of the renormalization factor in fc-space. Each renormalization factor calculated 
from the energy-range [-0.5,0.5] (eV) around the Fermi level except for that with index "wide" which is calculated from the 
energy-range [-10,10] (eV) around the Fermi level. 



z 


Average 


r 


X 


M 


R 


^wide 


Averaj 


;e 


r 


X 


M 


R 


xy 


0.55 


0.54 
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0.77 


zx 


0.72 




0.73 
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TABLE III: Transfer integral and its components for the tig orbitals of the V sites in the bulk SVO, t mn (R x ,Ry,R z ), where 
m and n denote symmetries of tig orbitals. Units are given in meV. t LDA is the expectation value of the KS-Hamiltonian 
for the Wannier function : t LDA = (0 L \H LDA \4> L ) (Eq. ©). t GWA is the static transfer integral in the GWA : t GWA = 
(<?i i |'H LDA + Z{— V xc + ReE)|<^ L ). " V" xc " is the expectation value of the exchange-correlation potential for the Wannier function 
. „yxc,,_ (0 L \V* C \(/) L }. "(Eh + A£ L )(0)" is the expectation value of the constrained self-energies at ui = for the Wannier 
function : "(Eh + AEx)(0)" = (<j> L \C£ H + AY, l )(uj = 0)\<j) L ). t H is the static transfer integral without double-counting : 
t H = {(f) L \n LDA + Z H (-V* C + ReH H )\4> L ) (Eq. t HA is the static transfer integral including the correction of the in- 

dependence of U : t HA = {4> L \H WA + Z HA (~V* C + Re(E H + AE L ))\cf> L ) (Eq. fl32)). Figure in the bracket in [0,0,0] is the 
value of the constant shift from the Fermi level of the high-energy space. 
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pair becomes about 50% larger than that of the LDA. In 
other words, the absolute value of the "effective" nearest 
transfer integral is about 130 meV lower than that of the 
"bare" nearest transfer integral given by £ LDA — V xc be- 



cause of the screening effect of the high-energy degrees 
of freedom. 

We now consider the self-energy effect arising from the 
w-dependence of the partially screened Coulomb inter- 
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TABLE IV: Effective Couiomb interaction between the two electrons for all the combinations of V tig orbitals in SVO (in 
eV). v and J v represent the bare Coulomb interaction/exchange interactions respectively. (7(0) and J(0) represent the static 
values of the effective Coulomb interaction/exchange interactions (at u = 0). The index 'n' and 'nn' represent the nearest V 
site [1,0,0] and the next-nearest V site [1,1,0] respectively. 



SVO 


V 


U(0) 


J V 


J(0) 




xy yz zx 


xy yz zx 


xy yz zx 


xy yz zx 


xy 
yz 
zx 


16.03 14.86 14.86 
14.86 16.03 14.86 
14.86 14.86 16.03 


3.51 2.52 2.52 

2.52 3.51 2.52 
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0.55 0.55 
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0.55 0.55 


0.47 0.47 
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FIG. 5: (Color online) Electronic band structure of SVO, 
where we calculate E H = —G H W for t2 g Wannier orbitals 
instead of V* c according to Eq. (|16[) . and subtract the aver- 
age of the self-energy from the low-energy bands. The zero 
energy corresponds to the Fermi level. 



action. Figure [6] is the w-dependence of the partially 
screened Coulomb interaction between the V ti g orbitals. 
As is the case with the constrained self-energy T, H , the 
w-dcpendence of the partially screened Coulomb interac- 
tion W p is weak for small u>(< 10 eV). We summarize 
the static limit of the partially screened Coulomb inter- 
action in Table IIV1 We renormalize the cj-dependence 
of the screened Coulomb interaction neglected in Table 
IIVI into the one-body part by following Eq. (f22j) . Fig- 
ure [7] shows the w-dcpcndcnce of the correlation term of 
AE L = ~G L (W - W u ). Because A£ L originates from 
the low-energy degrees of freedom, the overall slope in 
the w-dcpcndencc of AE L is basically similar to G L W il- 
lustrated in Figl^a). Small dip originating from L-L po- 
larization G L G L is observed near the energy eigenvalue 
of the low-energy band. 

Figure [5] illustrates the band structure where the w 
dependent part of the partially screened Coulomb inter- 
action is renormalized into the one-body part by follow- 



ing Eq. This band not only includes the effect of 

Z A = (1 - ^^y 1 but also includes the effect of Z H . 



In Table IIII[ the nearest transfer integral renormalized 
by the w-dependence of U, t HA , is substantially reduced 
from t H . This reduction of the transfer integral is mainly 
caused by the renormalization factor Z HA . 

It is remarkable that the resultant band structure of 
the effective model is nearly the same as that of the LDA. 
TableUshows that the final band width is 2.56 eV in com- 
parison to the LDA result of 2.58 eV. The two corrections 
originated from different effects, namely the increase in 
the band width arising from Y, H and the reduction aris- 
ing from AS L are roughly compensated and results in 
t HA similar to the LDA estimate. 



B. FeSe and FeTe 

FeSe^Tei-j; is the simplest iron-based 
superconductor—. It shows superconductivity with 
the transition temperature T c ~ 10 K at ambient 
pressur o 71 ' 72 and T c ~ 37 K under pressure (7 GPa) 7 -^. 
FeTe indicates an antiferromagnetic ordered moment 
~ 2.0-2.25^b at a Bragg point (7r/2, 7t/2) in the extended 
Brillouin zon o 72 ' 74 , whereas magnetism is not observed 
in FcSc. In this subsection, we derive the effective model 
of the two compounds and discuss how orbital-dependent 
effective parameters affect low-energy properties. 

Figure M (a) shows the band structure of FeSe in 
the DFT/LDA. As is the general case in iron-based 
superconductors 4 ^, ten states having strong Fe 3d charac- 
ter form a band near the Fermi level. The band is about 
4.5 eV in width, and is occupied by twelve electrons per 
unit cell. Small electron pockets are found around the M 
point and hole pockets are around the T point. Below the 
d bands, six states exist between —6 and —3 eV. They 
consist of mainly Se Ap orbitals. 

We regard the ten states having Fe 3<i character as 
the low-energy states, and construct the effective model. 
Figure [TU] shows the isosurface of the maximally localized 
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FIG. 6: (Color online) Frequency dependence of partially screened Coulomb interaction U of the ti g orbitals of the V sites. 
Here, v is the bare Coulomb interaction, and v is the average value of the bare Coulomb interaction of the V tig orbitals. 
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FIG. 7: (Color online) Correlation part of AE L of the ti g 
Bloch functions of SVO at T = (0, 0, 0) in the Kohn-Sham 
gauge. 



Wannier functions associated with the ten states. They 
are spatially extended because of hybridization between 
the Fe 3d atomic orbital and adjacent Se Ap atomic or- 
bitals. Since the strength of hybridization depends on 
the orbital, spread of the Wannier functions is orbital- 




FIG. 8: (Color online) Electronic band structure of SVO, 
where we calculate T, H + AE L for tig Wannier orbitals in- 
stead of V xc according to Eq. I|29p . and subtract the average 
of the self-energy from the low-energy bands. The zero energy 
corresponds to the Fermi level. 



dependent (Tabic W\ . This is in sharp contrast with 
the t2g-Wannier orbitals of SVO in the previous sec- 



tion. The yz/zx and 



y 



orbitals are delocalizcd 
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FIG. 9: (Color online) Electronic band structures of FeSe obtained by (a) the LDA and (b) the GWA and calculated according 
to (c) Eq. (|16l) (double-counting-less of the low-energy space) and (d) Eq. (|29[) (renormalizing the cj-dependence of W p ). The 
zero energy corresponds to the Fermi level. 



(a) dxy (b) dyz 




FIG. 10: (Color online) Isosurface of the maximally localized 
Wannier functions ±0.02 a.u. for the Fe 3d orbitals in FeSe. 

compared to the xy and 3z 2 — r 2 orbitals. This trend 
is observed not only in FeSe but also in FeTe and other 
iron-based superconductors^. (The xy axes in our con- 



vention are along the unit vectors of the cell containing 
two Fe atoms^. We abbreviate the 3d orbitals such as 
d xy as xy, unless confusions occur.) 

Orbital dependence is also observed in the partial den- 
sity of states (pDOS). Figure ITT1 presents the pDOS of 
FeSe resolved by the Wannier functions. The x 2 — y 2 or- 
bital has large density of states at the Fermi level in the 
LDA. A dip (pseudogap) is seen at about 0.3 eV above 
the Fermi level. Table fVTl shows the occupation number, 
where the 3z 2 — r 2 orbital has the largest value (nearly 
3/4- filling), while the x 2 —y 2 and yz/zx orbitals are about 
half-filling. We will discuss the results other than LDA 
later. 

Figures |H] (b) is the band structure of FeSe in the 
GWA, where we add the self-energy to the low-energy 
states after removing the LDA exchange-correlation po- 
tential. The band dispersion is similar to the LDA result, 
though the bandwidth becomes about 10% smaller by the 
self-energy effect. This trend is more clearly seen in the 
pDOS. In Fig. [IT] (b), we can confirm that the pDOS is 
similar to that of LDA except for overall band narrow- 
ing. The occupation number of each orbital is nearly the 
same as the LDA value. 
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FIG. 11: (Color online) Partial densities of FeSe resolved by the Wannier function of the Fe 3d orbital: (a) the LDA and (b) the 
GWA and (c) according to Eq. (116[) (double-counting-less of the low-energy space) and (d) according to Eq. ([29} (renormalizing 
the aj-dependence of W p ). The zero energy corresponds to the Fermi level. 



TABLE V: Spread of the Wannier functions (in A 2 ) defined 
by quadratic extent. 



Spread 


xy 


yz 


o 2 2 

Sz — r 


zx 


2 2 

x -y 


FeSe 
FeTe 


1.66 
2.35 


1.94 
2.76 


1.58 
2.01 


1.94 
2.76 


2.36 
2.43 



Now we turn to the results of effective models. Figure[9] 
(c) shows the band structure for Eq. (|20[) , where the self- 
energy correction from high-energy states is included. As 
is observed in SVO, the band is substantially wider than 
that of LDA. In contrast with SVO, however, the band 
dispersion is qualitatively different from the LDA band 
structure. Both the small electron pockets around the 
M point and hole pockets around the T point disappear. 
We note here that this band structure does not include 
the self-energy effect from the low-energy space, hence it 
should not be compared to experimental measurements. 



TABLE VI: Occupation number of the Wannier functions of 
FeSe, where the sum of the occupancy is 6. 



Occ. Num. 


xy 


yz 


o 2 2 

3z — r 


zx 


2 2 

x -y 


LDA 


1.23 


1.05 


1.57 


1.05 


1.11 


GWA 


1.22 


1.03 


1.53 


1.03 


1.19 


n H 


1.18 


1.30 


0.86 


1.30 


1.36 




1.19 


1.01 


1.43 


1.01 


1.35 



The partial density of states and occupation number cor- 
responding to Fig. [9] (c) are shown in Fig. QT] (c) and 
Table EH respectively. The psudogap observed in both 
LDA and GWA disappears. The x 2 — y 2 orbital has a 
weight down to ~ —6 eV below the Fermi level, while 
the xy component extends up to 3 eV, which results in 
the band widening. The occupation number of the yz, 
zx and x 2 — y 2 increases, whereas that of the 3z 2 — r 2 
decreases substantially to 0.86 from the LDA value of 
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FIG. 12: Real parts of correlation self-energy of the 3d Bloch functions at T — (0, 0, 0) in the Kohn-Sham gauge. 



25 
20 
15 
10 
9 5 
& 
-5 
-10 
-15 
-20 



Kohn-Sham gauge (at r) 



ImE (correlation part) 




4 
3 
2 
1 

I 
"""-1 
-2 
-3 
-4 



-30 -20 



-10 10 

co(eV) 



20 30 



ImE (correlation part) 




-30 -20 



-10 10 

co(eV) 



20 30 



FIG. 13: Imaginary parts of correlation self-energy of the 3d Bloch functions at T = (0, 0, 0) in the Kohn-Sham gauge. 



1.57. Consequently, the order of the occupation numbers 
is completely different from the LDA and GWA. 

From now on, we focus on the effective model param- 
eters obtained from H HA in Eq. (f3"2")l . because they are 
the best parameters which should be used in low-energy 
solvers. In "H HA ", the oj dependence of the screened 
Coulomb interaction is renormalized to the one-body 
part. The static value of the on-site screened Coulomb 
(U) and exchange (J) interactions in FeSe are presented 
in Table I Villi We can see that the U strongly depends 
on the orbital. Inclusion of its w-dependence using Eq. 
(pl2"j) yields the band structure shown in Fig. [5] (d). The 
band becomes narrower by the dynamical effect of U. 
The bandwidth is about 4 eV, which is smaller than the 
LDA bandwidth and close to the GW one. The pDOS is 
shown in Fig. [TT] (d). The yz/zx and 3z 2 — r 2 compo- 
nents have a peak near the Fermi energy. This explains 
why the numbers are quite different between LDA, GWA, 
H H and n HA . 

These results are understood from the transfer in- 
tegrals and on-site energies. The transfer integrals, 



t mn (R), are summarized in Table IVlTl In the tables and 
this subsection, the symmetry of the d orbitals is de- 
noted as the number; 1 for xy, 2 for yz, 3 for 3z 2 — r 2 , 
4 for zx, and 5 for x 2 — y 2 orbitals^. The on-site en- 
ergy listed in the column for (R x , Ry, R z ) — (0, 0, 0) and 
the nearest-neighbor transfer integral listed in the col- 
umn for (1/2,-1/2,0) arc extracted in Figs. IH1 and IT51 
respectively. 

It is useful to discuss the behavior of y. H to under- 
stand the parameter values for H HA . In comparison to 
the LDA and H HA approximation, the on-site energy of 
the x 2 — y 2 is lower in TL H , and the transfer integral for 
x 2 ~ V 2 (778 meV) is much larger than that in LDA (55 
meV). Both of them widen the pDOS of the x 2 — y 2 or- 
bital. The on-site energy of the 3z 2 — r 2 is higher in H H , 
which results in reduction of the occupation number. The 
next-nearest transfer integrals (R X7 R y , R z ) = (1,0,0) 
are comparable with the nearest-neighbor ones in the 
LDA. Especially, the (m, n) = (zx, zx) component, t'lf A 1 
is nearly the same as t\f A , which makes the system frus- 
trated. These strong tendencies in H. H arc substantially 
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TABLE VII: Transfer integral and its components for the 3d orbitals of the Fe sites in the FeSe, t mn (R x , R y , R z ), where 
m and n denote symmetries of 3d orbitals. Units are given in meV. f LDA is the expectation value of the KS-Hamiltonian 
for the Wannier function : LDA = ((l> L \'H Lr ' A \(f> L ) (Eq. J2j). is the static transfer integral without double-counting : 
t H = (0 i |'H LDA + Z H (-V XC + ReE H )|<?i L ) (Eq. (O). t HA is the static transfer integral including the correction of the in- 
dependence of U : t HA = (4> L \n LDA + Z HA (-V™ + Re(E" + AE L ))\<j> L ) (Eq. fl32)). Figure in the bracket in [0,0,0] is the 
value of the constant shift from the Fermi level of the high-energy space. 
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TABLE VIII: On-site effective Coulomb interaction between the two electrons for all the combinations of Fe 3d orbitals in 
FeSe (in eV). (7(0) and J(0) represent the static values of the on-site effective Coulomb interaction/exchange interactions (at 
u) = 0). 
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FIG. 14: On-site potential of the Wannier function of FeSe. 
Figure in the bracket is the value of the constant shift from 
the Fermi level of the high-energy space. 



> 0.6 

u 
cc 

1 0A 

0.2 

o 



I 

> 0.6 
CC 

S0.4 

LU 

0.2 





LDA 



(a)xy 



(0 3z 2 -r 2 



n 



FeSe If | [1/2,-1/2,0] 



LDA jji JJHA 



01 

> 0.6 
g 0.4 
0.2 

- 



(b)yz 



LDA 



s 0.8 

QJ 








> 0.6 
KD 








% °- 4 








0.2 




n 


n 


n 


i 


LDA 


(d)x 2 -y 2 





n 



FIG. 15: Nearest-neighbor transfer integral of the Wannier 
function of FeSe. 



weakened in H , but still a similar difference remains 
relative to the LDA result. 

In the H HA case, the order of the on-site energies is 
partially changed from the T-L H case. The crystal-field 
splitting is reduced and the transfer integrals change in 
both magnitude and order, which brings qualitative mod- 
ification of the band dispersion. Compared to the LDA, 
the on-site energy of the x 2 — y 2 is lower by 0.2 eV and the 
occupation number increases, while the 3z 2 — r 2 orbital 



is shifted upward by as large as 0.4 eV and the occupa- 
tion number decreases. As is the case with the LDA, the 
nearest-neighbor transfer integral for (m,n) = (xy,xy), 
t^ A is the largest in magnitude. The value is —460 meV 
and is about 10% larger than that of the LDA (—409 
meV). The (to, n) = [zx, zx) transfer, t'H A , is about 30% 
smaller than t^ A . The transfer integrals t HA except for 
the yz/ zx orbitals become larger than those of the LDA, 



16 



TABLE IX: Occupation number of Wannier functions of FeTe, 
where total occupancy is 6. 



Occ. Num. 


xy yz 


o 2 2 

3z — r 


zx 


2 2 

x -y 


LDA 


1.27 0.96 


1.58 


0.96 


1.23 


GWA 


1.25 0.99 


1.54 


0.99 


1.24 


U H 


1.28 1.52 


0.86 


1.52 


0.82 




1.21 1.25 


1.35 


1.25 


0.94 



and the increase is substantial for the 3z 2 — r 2 and x 2 ~y 2 . 
Since the x 2 — y 2 orbital is primarily responsible for the 
strong correlation effects with the Mott physic o 47 ' 48 , the 
increase in its transfer may cause some reduction of the 
correlation effects. 

We next show results for FeTe. Figure [TBI fa) shows the 
band structure of FeTe in the LDA. The ten-fold Fe 3d 
bands partially entangle with the Te bp valence bands, 
while the low-energy bands are similar to that of FeSe. 
In fact, the density of states of FeTe for the Fe 5d bands 
shown in FigurcfTTlfa) reveals that the x 2 — y 2 orbital has 
the largest density of states at the Fermi level followed by 
the yz/zx orbitals similarly to the case of FeSe discussed 
above. In case of FeTe, these two orbitals even have 
peaks near the Fermi level. In Table IIX1 we see that the 
occupation number is also similar to that of FeSe. 

Figures [TBI fb) andfTTKb) show the band structure and 
the density of states of FeTe in the GWA, respectively, 
where we disentangle the Fe 3c? Wannier bands from the 
whole KS-band structure 1 -^ and calculate the self-energy 
only for them. The distribution of the density of states 
and the occupation number are again nearly the same as 
those of the LDA and similar to those of FeSe (see Table 
EH). 

We now show in Figs. [16] (c) and (d) the band structure 
of FeTe after including the self-energy effect described by 
H H introduced in Eq. d20j) and H HA in Eq. ([32]), re- 
spectively, excluding the double counting of the electron 
correlation. As is the case with FeSe, the band structures 
are substantially different from that of the LDA. The cor- 
responding densities of states are also shown in Figs. [T7] 
(c) and (d). We summarize the transfer integrals of FcTc 
without the double counting of the electron correlation 
in Table H 

Now we focus on the effective model parameters ob- 
tained from M HA in Eq. ([32]) . The transfer integrals 
t HA except for the yz/zx orbitals become larger than 
those of the LDA, and the increase is substantial for the 
3z 2 — r 2 and x 2 — y 2 , similarly to FeSe. 

Geometrical frustration effects measured by the ratio 
between the next-nearest-ncighbor (i') to the nearest- 
neighbor (t) transfers also show a tendency similar to the 
case of FeSe: The frustration \t' /t\ is remarkably sup- 
pressed from 30/71 to 56/279 for the diagonal transfer 
between two x 2 — y 2 orbitals. Although this reduction 
looks dominating the magnetic stability, the frustration 
for the 3z 2 — r 2 orbital increases because t and t' both be- 
come more than twice of the LDA parameters and have 



similar amplitudes (t = —99 meV and t' = —154 eV). 
All of these corrections from the LDA results are qual- 
itatively similar to the case of FeSe. The increase of t' 
between the 3z 2 — r 2 and x 2 — y 2 orbitals to 66 meV is 
also noticeable, which is different from the case of FeSe. 

The stripe order is in general stabilized when t' be- 
comes comparable to t in amplitude, while the reduction 
of \t' /t\ may stabilize longer-period antiferromagnetic or- 
der such as up-up-down-down structure 1 ^ as is observed 
as the double stripe order in FeTe before the simple stag- 
gered G-type antiferromagnetic order becomes stabilized. 
The previous solution for the ab initio low-energy mod- 
els by the variational Monte Carlo method 4 ^ based on 
the LDA band structure has shown nearly degenerate 
ground states of the stripe and double stripe structure. 
By considering the overall reduction of the frustration 
particularly for the pair at the x 2 — y 2 orbitals, it is an 
interesting issue to study whether the double stripe phase 
becomes the unique ground state in the present ab initio 
model as in the experimental observation. 

The effective Coulomb interactions for FeTe obtained 
from the cRPA are listed in Table IXll The on-site U, the 
exchange J and the nearest-neighbor diagonal Coulomb 
interaction V are very similar to those obtained previ- 
ously^ by using 'H LDA (namely, the LDA band structure) 
and the differences are all within 10 % and mostly below 
0.2 eV. Similarly to the previous results 4 ^, the overall 
correlation amplitudes are slightly smaller than those of 
FeSe. 

A sharp contrast emerges in the comparison of the ab 
initio models for FeSe and FeTe: The contrast is found 
in orbital level shifts (or the on-site potential shifts) of 
the x 2 — y 2 orbital when we compare with the LDA band 
structure, as we see in the comparison between Figs. [Jj] 
and 1181 For the present FeTe ab initio model, although 
the relative shift of the x 2 — y 2 orbital level is not remark- 
able in Fig[l8] the occupation number of the x 2 — y or- 
bital is reduced from the LDA value 1.23 to 0.94 as we list 
in Table HXl On the contrary, FeSe shows a completely 
opposite behavior, where the occupation number of the 
x 2 — y 2 orbital increases from the LDA value 1.11 to 1.35, 
because of a substantial downward shift of the x 2 — y 2 
orbital level relative to other orbitals. 

We now elucidate the origin of this conspicuous down- 
ward shift of the x 2 — y 2 orbital of FeSe. We find that 
this downward shift is even stronger in t H (actually the 
x 2 — y 2 orbital becomes more than 0.5 eV lower than 
any other Fe 3d orbital as one sees in Table IVlTl In con- 
trast, the occupation number of the xy orbital, which 
has the band-insulator-likc pDOS, is almost unchanged 
from that of the LDA (~ 6/5 = 1.2). Wc then show the 
band structure, the density of states, and the transfer 
integral calculated without the exchange-correlation po- 
tential, t° = (<j) L \H LDA - V xc \cj) L ), in Figs. [JJ and HO] 
and Table IXHi respectively. Overall behavior is nearly 
the same as that of "t H " except for the uniform reduction 
of the band width. The total width of the Fe 3d bands 
in is about 80% larger than that in "£ H " . This is 
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TABLE X: Transfer integral and its components for the 3d orbitals of the Fe sites in FeTe, t mn (R x ,R y , R z ), where m and 
n denote symmetries of 3d orbitals. Units are given in meV. t LDA is the expectation value of the KS-Hamiltonian for the 
Wannier function : LDA = (<j> L \H WA \4> L ) (Eq. ©). t H is the static transfer integral without double-counting : t H = 
(0 L |"H LDA + Z H (-V* C + ReE H )|0 L > (Eq. pCi]l). t HA is the static transfer integral including the correction of the ^-dependence 
of U : t HA = {<f} L \n WA + Z HA (-V* C + Re(E ff + AE 1 '))^ 1 '} (Eq. (E2J)). Figure in the bracket in [0,0,0] is the value of the 
constant shift from the Fermi level of the high-energy space. 
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FIG. 16: (Color online) Electronic band structures of FeTe obtained from (a) LDA, (b) GWA, (c) Eq. (H H , excluding 

double counting in low-energy space) and (d) Eq. ()32|) (H HA , renormalizing the cj-dependence of W p in addition to H H ). The 
zero energy corresponds to the Fermi level. 



because the self-energy £ only weakly depends on the 
wave number and orbital in contrast to AY> L . Therefore 
T, H uniformly shrinks the band width from "t°" . 

Since the downward shift of the x 2 — y 2 orbital level 
in t° is even more enhanced than t H , we find that the 
essence of the shift is contained in the procedure of sub- 
tracting the exchange correlation potential of the LDA: 
When the effect of the exchange correlation would be re- 
moved, the x 2 — y 2 orbital would be located at a level 
much lower than other orbitals, because its stronger hy- 
bridization with the chalcogen p orbitals contributes to 
the exchange correlation potentials. By contrast, in FeTe, 
because the yz/zx orbitals have the strong hybridization 
with the chalcogen p orbitals (see Table |Vj , the reduc- 
tion of the on-site potential of the yz/zx orbitals is the 
largest. Such weakly localized orbital, whose effective on- 
site interaction is relatively weak, has weaker self-energy 
compared to the other 3d orbitals. 

More strictly speaking, the contribution of the d elec- 
trons to the exchange correlation potential is more or less 
proportional to U in amplitude and the sign is negative. 
Therefore, if one removes this contribution, the level is 
pushed up more or less proportional to U . Since x 2 — y 2 



orbital has the weakest U for FeSe, the upward level shift 
is the smallest for x 2 — y 2 , meaning that it moves down- 
ward relative to other orbitals. This interpretation is in 
accord with the interpretation above by the largest hy- 
bridization of the x 2 — y 2 orbital, because the largest hy- 
bridization makes the largest Wannicr spread and hence 
the weakest U. The level shift is indeed in the order of U 
both for FeSe and FeTe and in the order of the Wannicr 
spread as well. Then the level shift of xy, yz/zx, 3z 2 — r 2 
and x 2 — y 2 orbitals of FeSe from the LDA result to the 
result by H HA are - +0.3, -0.05, +0.4 and -0.2 eV, 
respectively, as we see in Table I VIII Therefore, when we 
solve by using the present effective low-energy model with 
the single-particle part H HA instead of H LDA , the orbital 
levels should be effectively shifted with these amounts. 

Since the correlation effect is governed by the most 
correlated x 2 — y 2 orbital while the effective screened 
Coulomb interactions are nearly the same between the 
previous^ and the present models, such a downward shift 
of the x 2 — y 2 orbital may increase the filling of the x 2 — y 2 
orbital, and cause the depinning from half filling. This 
may destroy the antifcrromagnctic order as in the exper- 
imental observation. 




FIG. 17: (Color online) Partial densities of states for FeTe resolved by Wannier functions of Fe 3d orbitals calculated from (a) 
LDA, (b) GWA, (c) Eq. ([16} (H H , excluding double counting in low-energy space) and (d) Eq. ([29)l ("H HA , renormalizing the 
cj-dependence of W p in addition to 1-L H ). The zero energy corresponds to the Fermi level. 
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FIG. 18: On-site potential of Wannier orbitals for FeTe. Fig- 
ure in the bracket is the value of the constant shift from the 
Fermi level of the high-energy space. 



Here, we note again that the band structure and the 
density of states shown in Figs. [51 [TT1 [TBI and ITH should 
not be taken as the properties that can be directly com- 
pared with the experiments, because the experimentally 
accessible quantities such as the Fermi surface structure 
and the spectral weight are obtained only after solving 
the effective low-energy models. 



IV. SUMMARY 

We have proposed an improved scheme for ab initio 
derivation of the low-energy effective models. Our for- 
malism is free from the double counting of Hartree and 
Fock contributions from the low-energy space, and the 
LDA exchange correlation is replaced by the GW self- 
energy. Moreover, the derived effective model is reduced 
to a static one, where the w-dependence of the screened 
Coulomb interaction is taken into account and rcnormal- 
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TABLE XI: Bare and effective Coulomb interactions between two electrons for all the combinations of Fe 3d orbitals in FeTe 
(in eV) . Here, v and J v represent the bare on-site and exchange Coulomb interactions, respectively. The static limit (u —> 0) of 
the effective on-site and exchange Coulomb interactions are denoted by U(0) and J(0), while v n and V(0) represent the bare 
and effective nearest-neighbor Coulomb interactions, respectively. 
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FIG. 19: (Color online) Electronic band structures calculated without exchange-correlation potential of FeSe and FeTe. The 
zero energy corresponds to the Fermi level. 



ized to the one-body part as a self-energy. We have ap- 
plied this formalism to transition-metal oxide SVO, as 
well as to iron-based superconductors FeSe and FeTe. We 
have found there are two opposite effects, namely the in- 
crease in the band width arising from Y, H and the reduc- 
tion arising from AS L . In SVO, these effects arc roughly 
compensated, and remarkably, the resultant band width 



and the dispersions are nearly the same as that of the 
LDA (the band width of LDA: 2.58 cV, H HA : 2.56 eV). 
On the other hand, in the non-degenerate multi-band 
systems such as FeSe and FeTe, though the band widths 
are also similar to the LDA results, the momentum and 
orbital dependent self-energy effects yield modifications 
of the resultant band structures. For the effective low 
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FIG. 20: (Color online) Partial densities of states calculated without exchange-correlation potential for FeSe and FeTe resolved 
by Wannier functions of Fe 3d orbitals. 



energy model for FeSe given by H , the on-site po- 
tential of the x 2 — y 2 orbitals, which has the strongest 
hybridization with the Se 4p orbitals and thus the weak- 
est on-site interaction among the Fe 3d orbitals, is sub- 
stantially lowered from that of the LDA. This may make 
the occupation number of the x 2 — y 2 orbital away from 
half filling and destroy the antiferromagnetic order as in 
the experimental observation. In contrast, the lowering 
of the on-site potentials is the largest in the yz/zx or- 
bitals in FeTe, and the resultant occupation number of 
the x 2 — y 2 orbital remains close to half filling. 

By our formalism, two major drawbacks in the deriva- 
tion of the low-energy effective model in the literature are 
removed. The effects of Y, H and AS L compensate and 
the resultant band width is, rather accidentally, nearly 
the same as that of the LDA, while the structure of the 
band and the on-site potential of the nonequivalent or- 
bital quantitatively change in some systems such as FeSe 
and FeTe. 

Even with this improvement, we still have future is- 
sues to be resolved: In the practical calculation of the 
low-energy effective model by low-energy solvers, one of- 
ten neglects the farther-neighbor effective interaction ow- 
ing to heavy costs of the numerical computation. This 
long-range part makes the effect of the on-site interaction 
weaker. The polaronic screening effect from the lattice 
also weakens the effective Coulomb interaction between 
the electrons. When simplified models are used in the 
low-energy solvers, these correction may have small con- 
tributions. Counting of these effects to the downfolding 
formalism are left for future problems. 

Another open issue is to consider possible refined esti- 
mates of the GW self-energy arising from the high-energy 



contributions. A self-consistent GW scheme instead of 
the one-shot self-energy has been employed by Kutepov 
et ali^ with the combination of the solver based on the 
dynamical mean field theory. In the full GW calcula- 
tions, it is known for a long time that the self-consistent 
GW calculations give worse agreements with the exper- 
imental estimates of the gap amplitudes than those of 
the one-shot GW. This may be attributed to the ver- 
tex corrections ignored in the GW approximation that 
roughly cancels the difference between the one-shot and 
self-consistent estimates. Counting both of the vertex 
corrections and self-consistency is, though expected to 
cause minor differences in this constrained self-energy, 
left for future studies. 
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TABLE XII: Transfer integral between Fe 3d Wannier orbitals for FeSe and FeTe, t mn (R x , R y , R z ), where m and n denote 
symmetries of 3d orbitals. Here, Fock terms are not included. Units are given in meV. Transfer integral calculated without 
exchange-correlation potential is given as t° = (4> L \'H LUA — V xc \<j) L ). Figure in the bracket in [0,0,0] is the value of the constant 
shift from the Fermi level of the high-energy space. 
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